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Abstract 

The stopping power and energy loss rate of charged particles traversing a two-dimensional Dirac 
plasma is investigated. The Dirac plasma considered here models a solid state system, recently 
realized graphene monolayer, where the conduction electrons obey the Dirac-like equation and 
exhibit a linear in momentum dispersion relation. Theoretical work presented here is based on the 
the dielectric response function and the dynamical structure function within the Random-Phase- 
Approximation (RPA). 
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I. INTRODUCTION 

Recent fabrication of single layers of graphite (graphene monolayer) {l| has posed the inter- 
esting question of charged particle energy loss in solid media where the conduction electrons 
behave distinctly different from those found in ordinary semiconductor heterostructures. 
In a two-dimensional electron gas (2DEG) realized in semiconductor heterostructures, the 
conduction electrons behave as ordinary electrons with a parabolic dispersion relation. In 
contrast, the dispersion relation obeyed by conduction electrons in graphene monolayer is 
linear in momentum. This occurs due to the unique crystal structure of graphene which is 
a two-dimensional (2D) honeycomb lattice of carbon atoms. Quantum-mechanical hopping 
between the sublattices of graphene leads to the formation of two energy bands, and their 
intersection near the edges of the Brillioun zone yields the conical energy spectrum. As 
a consequence, the dispersion relation of electrons and holes bands is linear near K' 
points of the Brillioun zone which is given by = hvpk. Hence, the conduction elec- 
trons, known as Dirac electrons, behave as massless particles with the effective speed of 
light vp — lO^m/s In a graphene monolayer, electron transport is essentially 

governed by Dirac's (relativistic) equation, rather than the usual Schrodinger equation for 
nonrelativistic quantum particles. The relativistic behavior of graphene was first predicted 
by P. R. Wallace [3]. It thus provides a unique opportunity to study relativistic quantum 
dynamics in condensed-matter systems. The Dirac-like gapless energy spectrum was also 
confirmed recently by cyclotron resonance measurements in graphene monolayer jj]. 

In this work, we address energy loss experienced by charged particles as they traverse a 
non-local dynamic quantum plasma of Dirac electrons, electrons that obey the Dirac equa- 
tion. The energy loss calculation performed here requires that we determine the dynamic, 
non-local dielectric response function for the plasma of Dirac electrons. We determine the 
dielectric response function within the Random-Phase-Approximation(RPA) and employ it 
to calculate the dynamical structure function which yields the stopping power. 

Since the work of Nozieres and Pines 51], conduction electrons in solids have been treated 
as a solid state plasma, a gas of electrons in a neutralizing positive background. Moreover, 
electron energy loss spectroscopy has been an important probe of the dielectric response 
properties of solid state systems^]. This was brought into clear focus by the pioneering 
work of Ritchie j?] where the energy loss of a charged particle traversing a bounded solid 
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state plasma was considered. Dynamical processes in a solid state medium have been inves- 
tigated with a variety of spectroscopic tools which employ an almost monochromatic beam of 
electrons, photons, neutral atoms or ions which scatter inelastically from the medium under 
study Information concerning the electronic properties of the medium is obtained as a 
result of the interaction between the probes and the elementary excitations of the medium. 
Experimentally, this information is obtained by electron energy loss spectroscopy (EELS) 
and inelastic light scattering studies. Recently, energy loss spectroscopy of free-standing 
graphene films was performed 

The present paper is arranged as follows. In section II we give the formulation of the 
problem. The energy loss rate and stopping power of a 2D Dirac plasma is outlined. This 
theoretical investigation is based on the evaluation of dynamical structure function and 
energy loss function. In section III, we present the results and discussion. Concluding 
remarks are made in section IV. 



II. FORMULATION 



The system under consideration is a 2D quantum plasma of Dirac electrons i.e. Dirac 
electrons embedded in a uniform and rigid neutralizing background of positive charges. We 
begin by considering the low-energy electronic eigenstates for this system, given by 

i/oF(r)=£F(r). (1) 



The Dirac-like Hamiltonian in eq(l) is lOj 



Q ]^ — i]^ 

Ho = Vf[^ . ^ ^ \ = VFicr-xK + (^yky) (2) 

~\~ iky 

where ax, Uy are Pauli matrices and fi? is the two-dimensional Fermi velocity of Dirac elec- 
trons (with the characteristic velocity vp — 10^m/s)[10]. Conduction electrons in a Dirac 
plasma behave as massless Dirac particles with a linear dispersion relation egk = skvp [l| 
where s = ±1 indicate the conduction (+1) and valence (—1) bands, respectively, k is the 
two-dimensional wave vector and is given by k = k^i + kyj, here = k cos 6^, ky = k sin 6k- 
Hence, electron dynamics is modelled by the following Dirac equation 

vpicTxkx + ayky)Fsk{r) = eFsk{r). (3) 
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From here onwards we take h = c = 1 throughout this calculation. The wavefunctions 



appearing in eq.([3]) are given by 



io| 



Fs,k(r) =^F,,kexp(ik.r), (4) 



Fs',k+q{r) = ^Fs',k+q exp(z(k + q).r) (5) 



1 /e-^^'=' 

F.. = ^ „ ) (6) 



where is the area of the plasma sheet and 

1 

V2\ s 
and 

The wavevector and frequency dependent longitudinal dielectric response function within 
the Random-Phase-Approximation(RPA) can be expressed asj?! 

e(q,a;) = 1 - Ucn{q,uj) (8) 

where z/^ = 27re^//€g is the Fourier transform of two dimensional coulomb interaction, k is 
the background dielectric constant and 7r(q, cj) is the two-dimensional(2D) polarizability. 
The plasmon modes at finite wave vectors are given by the zeroes of the dielectric response 
function given by eq(IH])- Since, in this work, we are primarily interested in the energy loss 
due to plasmons, we consider the dielectric response function in the high frequency, long 
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wavelength limit (g — > 0) , which is given as 11 1 



Here Up = {gsgvG'^Ep/2K)^^'^ ^^q is the plasma frequency with Ep the Fermi energy. Equiv- 
alently, in terms of dimensionless variables, we have 

2'Krs fgsgvney/"^ X f^_J^\ ^^Q^ 



In the above expression, we have used the dimensionless variables x = q/kp, v = uj/Ep,.rs is 
the Wigner-Seitz radius, k is the background dielectric constant and gs = g^ = 2 being the 



spin and valley degeneracies respectively UJ]. Energy loss function is the basic parameter 
which accounts for the energy lost by the incident charged particle. It is also the parameter 
which is of central importance in energy loss spectroscopy (EELS) experiments. Energy loss 



function is defined as tlie imaginary part of the inverse dielectric response function. From 
the dielectric response function given by eq([9]) and introducing the infinitesimally small 
parameter t] we obtain 

1 iuj + ir]Y A B 



e{q,uj) 



(id + iT]y - cjM 1 - 



: + 

+ ir] u + Up./l 



(11) 

Evaluating the coefficients A, B and using Dirac's prescription lim^ >o = | -F ^7r5(,^), 

in eg (11 II) and then expressing the variables in terms of dimensionless parameters x v and 
we obtain the energy loss function for graphene monolayer as 



Im 



e(x, v) 



1/2 
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The interaction of charged particles with condensed matter system can be studied by means 
of the system's stopping power. The energy loss per unit path length is the stopping power 
S = —dE/dx and it accounts for the energy lost by an external charged projectile as it 
passes through matter. In the quantum mechanical framework, we consider an in-plane 
probe of charge ze, mass m and velocity fp, interacting with the many particle system under 
consideration by treating the incident particle state as a plane wave state. We calculate the 
probability that the point-like projectile loses an energy u; in a time interval dt in interaction 
with the 2D Dirac plasma. During the interaction impulse transfer is q. Here uo and q satisfy 
uj = q.Vp + q^ /2m. For a heavy incident particle, we neglect the recoil energy q^ /2m to 
obtain uj = q.v^.The stopping power is given by summing over energy difference weighted 
by transition rate W times the inverse projectile velocity, mathematically 



dE/dx = -l/vp^{ek - -q)W^. 
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We can expand energy loss in terms of Au = q^ /2m to obtain 12 1 



dx \dx J Q \dx J ^ 



(14) 



oo 

Here we also make use of the sum rule — J ujlm{l/e(q,uj)} dcu = 7rujp/2. To the lowest 



order, the stopping power for a 2D Dirac plasma given by eq (|T3|) can be expressed in terms 



of the dynamical structure functions S{q,Lj) as 

fdE\ 1 r rr, _,2^ fny^^V 



whereas the first order stopping power is given by 



The dynamical structure function is the auto-correlation function of the Fourier components 
of the particle density and accounts for longitudinal charge oscillations of the electron density 
relative to the positive background. From the polarizability, we can obtain the dynamical 
structure factor at finite temperature for a 2D Dirac plasma as 

S{(l,uj) = 7 -^(1 - exp(-/3t^))"^ [5{uJo - ujp) - 5{uJo + ujp)] (17) 

2g,g,e^neEF (l - ^) 



where cuq = cj/y 1 — oo'^/AEp, (3 = l/fc^T, fc^ is the Boltzmann constant, T is the absolute 
temperature, Ep = vpkp, kp = {A-Krie/ gsQvY^'^ is the Fermi wavevector and Ue is the electron 
density. Dynamical structure factor provides direct physical information about longitudinal 
excitations in the system. It is the measure of the density-density correlations of the system. 
Calculations of dynamical structure function reveal that the electromagnetic response of a 
Dirac plasma is substantially different from a 2DEG system which is essentially due to the 
distinctly different dispersion relation of Dirac electrons in a 2D Dirac plasma as compared 
to ordinary electrons in a 2DEG system. 

For a solid state Dirac plasma under consideration, we can classify the rate of energy 
loss and stopping power into two categories, one due to the particle-hole excitations and the 
other due to the plasma oscillations. So the total stopping power of the system is 

dE dEe—h dEpi . . 

dx dx dx 

Employing the following relation in eq IfTSll 

6{e -En-Vn cos(— a;o)) = 7^ 

^ .e{\Vn\-\e-en\) (19) 



^(^0-7^ arccos(^---^)) + 6{xo + arccos( ^ t/" )) 

/Vr Vri 2,71 V„ 
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For a single charged particle, z = 1, the lowest order stopping power, due to plasma oscil- 
lations, of a 2D degenerate Dirac plasma is given by 

.^^^-^ r . , ' ^ln(.+ v/^)arccos%).. 

(20) 



with y = qvp/upjl- 7^ 
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III. RESULTS AND DISCUSSION 



The dynamic and static response properties of an electron system are all embodied in the 
structure of the dielectric response function. We have employed the Random-Phase- Approx- 
imation (RPA) based dielectric response function to determine the structure function and 
the energy loss function of a 2D Dirac plasma. The two types of excitations in the medium 
responsible for energy loss are the collective and single particle excitations. Collective ex- 
citations (plasmons) occur at a higher energy compared to single particle (electron-hole) 
excitations. At a frequency ujq equal to Up, collective excitations (plasma oscillations) dom- 
inate for q'^Vp < Up and these exist as long as the wavevector remains less then a critical 
value qc- However, above this critical value plasmons undergo Landau damping generating 
electron-hole pair excitations. Therefore, energy loss rate and stopping power are affected 
by these two types of contributions, one due to the collective oscillations i.e. plasmons and 
the other due to the electron-hole excitations. The energy loss function taking into account 
both these mechanisms of energy loss is presented here. Collective modes with the onset of 
Landau damping change into electron-hole pairs. We have plotted the loss function given 
by eqf |T2l) due to collective excitations in Fig.(l). Here we essentially restrict the energy 
loss to plasmon excitations by plotting the result at a; = 0.1, 0.5, 0.01, 0.05, 0.001. The 
parameters used are: electron density Ue = 3.16 x 10^^/m^, vp = lO^m/ sec, k = 2.5(using 
Si02 as the substrate material), = 0.5. Note that the critical value of the wave vector qc 
and the corresponding critical value of x is 0.7 where the onset of Landau damping of plas- 
mons occur for the parameters considered here. Dynamical structure function shows delta 
function peaks at the plasma frequency. The peaks in the figure are the manifestation of 
existence of plasmons in the system. In Fig. (2), in order to include the contribution of single 
particle excitations, we have plotted energy loss function against u = uj/Ep for following 
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values ol X = q/kp — 1-0, 1.5, 2.0. The peaks in the energy loss function seen in the figure 
can be interpreted as occurring due to single particle excitations of the system. In Fig. (3), 
we have plotted energy loss function versus x = q/kp for various values of = LojEp. The 
relationship between energy loss rate dEpi/dt due to plasmons and velocity of the incident 
particle Vp is shown graphically in figure(4). We have plotted the energy loss rate versus 
the dimensionless velocity Vp/vp- Stopping power is shown in Fig. (5) where the energy lost 
by the incident particle is due to plasma oscillations. The stopping power decreases with 
the increase of the incident particle velocity. Furthermore, the lowest order energy loss rate 
and stopping power of the medium are found to be independent of the mass of the incident 
particle. 



IV. CONCLUSIONS 



An analysis of energy loss suffered by fast particles in their interaction with matter 
requires the calculation of energy loss rate and stopping power. In this work, we present a 
theoretical investigation of energy loss through the calculation of both the energy loss rate 
and the stopping power of a 2D solid state Dirac plasma. The motivation of this work is 
the recent realization of single layer of graphite, graphene monolayer. The quasiparticles 
in graphene monolayer are found to obey the Dirac equation with a linear in momentum 
dispersion relation. The work presented here is based on the the dielectric response function 
and the dynamical structure function within the Random-Phase- Approximation (RPA). The 
energy loss function taking into account the collective as well as single particle mechanisms 
of energy loss is presented here. In the collective excitation regime, energy loss function 
peaks at the plasma frequency, such peaks when detected by experiments, can also be used 
to identify energy at which plasmons occur in the medium. Furthermore, stopping power 
as a function of the incident particle velocity is also determined. Stopping power is found 
to decrease as the velocity of the incident particle is increased. Moreover, the lowest order 
energy loss rate and stopping power of the Dirac plasma are found to be independent of the 
mass of the incident particle. 
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